On Effective Electron Mass of Silicon MOSFET at Low Electron Density. 
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' The trial wave function method developed in Ref.sBII for the case of narrow s-band in a perfect 

, crystal is adapted for calculation of the density dependence of the effective mass and the Lande 

factor in a dilute two-dimensional electron system. We find that the effective mass has a tendency 
£SJ \ to diverge at a certain critical concentration, whereas the g factor remains finite. 
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PACS numbers: 71.30.+h, 73.40. Q\ 



As the temperature decreases, a dilute electron gas of highly mobile Si-MOSFET exhibits a strong drope in the 
resistance if the electron density n c is higher than a certain critical one n s > n c , and an inrease in resistance when 
£NJ I n s < n c a. In the vicinity of n c the resistance possesses scaling properties as a function of temperature and electron 
density. This instance from the very first caused researchers to consider the observed transition as a disorder-controlled 
1 quantum phase metal-insulator transition (MIT) and gave rise to a tide of similar investigations of other systems in 
which any change in the sign of derivative ^(n s ) was taken as an evidence of the occurrence of a quantum phase 

^: MIT - □ 

In the recent, experimental studyu of the screening properties of a two-dimensional electron system as a function 
■ of temperatureo, the interpretation was used in the terms of Ref.H , and a strong increase of the effective mass was 
found in Si-MOSFET as the electron density approached a value of ~ 0.8 * 10 11 cm -2 , which is close to n c in the best 
of investigated samples. A similar behavior of the cyclotron mass was observed in independent!!! experiments on the 
measurement of the temperature dependence of Shubnikov-de Haas oscillations. An analysis of the experimental data 
similar to that made inQ but performed in the opposite limit in the ratio of valley-splitting energy to temperature 
with the use of data of other experimental groups and samples from other sourcesffil, confirmed the versatility of the 
'"^J , m*{n s ) curve. 

The conclusion that should be made from recent experimental data is that the quantum phase transition observed 
in the most perfect MOSFET is rather the property of a pure disorder-free two-dimensional system. A— Qualitative 
i^i - theory of two-dimensional electron Fermi liquid in a state close to crystallization was presented in Ref.s&irE 2 ! Below 
a quantitative description of a two-dimensional paramagnetic electron liquid is proposed, adapting the trial wave 
function approach developed in Ref.saafor the case of a narrow s band in a perfect crystal. 

Let us assume that the ground state of an electron system with strong interaction in a regime close to crystallization 
can be described as an crystal with a great number of charge -carrying mobile defects. The real two-dimensional 
system will be replaced by a grid of lattice sites with a density of n s . An electronic wave function of the Wannier 

1 /2 

type 0[(r — g)n s ] , where the vector g specifies the position of a lattice site, will be associated with each site. The 
corresponding creation operator is a g . If each site were occupied by only one electron, the system would represent 
a perfect electron crystal. In fact, there is a certain probability depending on n s that in the ground state a site can 
be occupied by two electrons with opposite spins. The number of such states eventually determines the number of 
mobile excitations and, hence, the transport properties of the system. 
Bloch wave functions are constructed on the base of lattice sites: 



V>k(r) =n s 1/2 ^ea;p(ikg)0(r-g), (1) 



a k = ? V 1/2 I>^kg)at. (2) 
g 

The Hamiltonian of the system contains the electron kinetic energy and the electron interaction at one site: 

H = £k(4T a kT + a ki a ki) + 1/2 a l] a ll a sl a sV ( 3 ) 

k eon s g 

Here, £k = and £q is the static dielectric constant. We introduced a coefficient a into the interaction energy 

determined by the exact form of the wave function on the site and completely neglected the electron interaction on 
the neighboring sites. Processing to the limit of the gas of noninteracting electrons requires the modulation in Eg.(l) 
disappear and the coefficient a be a slowly varying function of n s , vanishing at nj 1 — > 0. We will neglect this weak 
dependence in the region of a low electron densities. 
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In Ref. 0, it was proposed that a many-body trial function be used for a ground state in the form 

v>=E^IRtIRi $ °> ( 4 ) 



or 



where G and T are the sets of sites occupied by electrons with spin up and down, repectivelly; and $0 is a vacuum 
state. It is convenient to express the function ip through operators of creation and annihilatin of Bloch waves and 
to take into account electron correlation by decreasing coefficient Aar in (||) by a factor of rf if the corresponding 
product implies the occurence of doubly occupied sites whose fraction equals v (0 < r\ < 1). The relation between v 
and 77 in the ground state was obtained iria. For our case, 



The probability that a single-particle state with a wave vector k is occupied undergoes a jump at k = kp by the value 

q=lSv(±-v). (6) 

Thus, the trial wave function describes a mixture of functions that corresponds to a fully occupied band (solid spin- 
ordered phase) and a paramagnetic electron liquid. The transition to the solid phase is continuous and is characterized 
by the parameter q (0 < q < 1): q = 1 in the paramagnetic electron liquid with weak interaction, and q = in the 
electron crystal. 

The mean value of the Hamiltonian given by (|^) in the state with given v equals 

(H) v = \n a qe F + 1/— n 3/2 , (7) 
2 £0 

where Ep is the Fermi energy of an equivalent number of electrons in the absence of interaction. According tell, the 
expression in Eq. (^) is minimized with respect to v with regard to Eq.(|^). A minimun of the Hamiltonian is attained 
at 



^ = i( 1 - ) V W = U — ^ ); 8 

which, according to (0), corresponds to 



5* 
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q- 1 = — = (9) 
m n s — n c i 



(10) 



Here, m* is the renormalizcd effective mass. In the same way, followinga, the Lande factor can be found as 

1 (n cl \l/2 



1 _ 1+ 2(^) 

U S ( 1 + ( ^)l/2)^j 



The simplest of way of generalization to the case of two valleys is in considering two parallel sublattices. In each of 
them the number of electrons equals n s /2 and the characteristic cell size is deminished compared to the single-valley 
case by a factor of (3. The coefficient (3 is determined by the ratio of Coulomb energies of inter and intravalley 
interactions. In the limit of two sublattices in one plane (3 = \/2 . In the case of two valleys, n c \ in Eq. (9,10) should 
be replaced by n C 2 = 2(3 2 n cl . 

A comparison of the curves obtained in this way with experimental results is shown in Figs. 1,2. Single fitting 
parameter n C 2 = 0.78* 10 11 cm ~ 2 has been used, which corresponds to a — 0.15. It is evident from the figures that the 
behavior of both the effective mass and the g- factor is reasonably described within framework of the proposed model, 
though the coefficient a is approximately twice as large as the value expected according to numerical calculationstj. 

It should be specially noted that the above considerations give no way of judging the spin state of the solid phase, 
because it is determined by the exchange interaction of electrons on neighboring sites. Moreover, in the immediate 
vicinity of the transition point, in the region where ((H(n s )) v — (H (n c )) ^nj 1 turns out to be smaller than the exchange 
energy of electrons on neighboring sites, the proposed description does not work in the paramagnetic electron liquid 
as well. Thus, the issues of the phase diagram in the immediate vicinity of the transition point and those of the spin 
structure of the solid phase remain out of scope of this consideration. In general, the approximation used is poorly 
controlled, and the rather good description of experiment still remains its only justification. 
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It is well known that the concentration n c that corresponds to a change of the sign of the derivative dR/dT strongly 
varies from sample to sample, depending on the disorder in the electron system under study. An impression is gained 
from experimental data that the transition point n C 2 measured-far different samples also somewhat varies. This fact 
can also be related to the effect of disorder, as was discussed inliillij. In strongly disordered electron systems n c ^> n C 2 
and the effect considered above, namely, the dramatic increase in the effective mass is not observed. In the most 
perfect of the electron systems studied, n c rs n C 2- 
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FIG. 1. Effective mass as a function of electron density. The solid line corresponds to Eq. |0| with n C 2 
Squares and circles correspond to the experimental data from Ref.u, and Ref.a, respectively. 
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FIG. 2. Effective g- factor as a function of electron density. The designations of experimental points are the same as in 
Fig. §. 
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